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ABSTRACT 


PEAPS (Passive Environmental ASW Prediction System) is 
a relatively unsophisticated model which accepts input source 
and receiver parameters and then predicts sound propagation 
characteristics. in an ocean environment, the corresponding 
transmission loss, and the probability of detection. The 
program was written for a programmable desk-top calculator 
for immediate deployment and operational testing aboard 
small ASW platforms. The program is also available in a 


form suitable for larger computers. 
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equation. 
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DEFINITION 
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INTRODUCTION 


A. BACKGROUND 

Shore-based prediction systems rely on complex acoustic 
models which require large computers and relatively long 
computation times. Their predictions are based on histori- 
cal meteorologic and oceanographic data for regions in which 
environmental conditions are assumed to be homogeneous. 
These systems use a single level of sonar performance as the 
basis for all range predictions. The value of such systems 
lies in their use as research tools and for pre-deployment 
planning: they can be considered to be strategic systems. 

With the advent of the programmable desk-top caiculator 
comes the possibility that a simple acoustic prediction model 
can be placed on board anti-submarine warfare (ASW) platforms. 
This makes it possible for the ASW platofrm to take advantage 
ШЕ Оп-5сепе bathythermographic (BT) data to obtain current 
range predictions based on the output of the platform's own 
sensors. In addition, an on-board system enables the ASW 
platform to obtain predictions under all electronic emission 
control (EMCON) conditions. The programmable desk-top cal- 
E E make possible environmental predictions at the 


tactical level, 
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В. PEAPS 
U Descr»ptren 

PEAPS (Passive Environmental ASW Prediction System) 
was designed for use in any programmable desk-top calculator 
or general-purpose computer with sufficient capacity for 
the program. 

The system is composed of four sub-programs: а 
speed of sound profile (SSP) generator, a ray trace generator, 
a transmission loss (TL) generator, and a probability of 
detection (PD) generator. The SSP and ray trace generators 
are combined as an SSP/Ray Trace program; the TL and PD 
generators are combined as a TL/PD program. These sub- 
programs are described in detail in subsequent chapters. 

Zee Assumptions and Constraints 

The system was written originally in a machine 
language for use on a small capacity, programmable, desk- 
Bop caiculator (the Hewlett-Packard Model 9810A with expanded 
Storage). The hardware used to develop PEAPS is described in 
Appendix A. Total program capacity of the calculator was 
2036 program steps with 100 data-storage registers. This 
calculator was used because it vas reletively inexpensive 
and readily accessible; it was reasoned that if a simple 
working model could be designed on a small-capacity, 
programmable calculator, more sophisticated models would be 
possible on larger programmable calculators. 

PEAPS is a relative unsophisticated system which 


uses approximation techniques to compute transmission losses 





under the assumptions of a flat-bottomed ocean and phase 
incoherent interference. Lateral variations in vertical 
temperature profiles are neglected. Speed of sound profiles 
are approximated by straight-line segments, with storage 
available for fourteen segments (fifteen input points) and 
five smoothed gradients. When available to the user deep 
ocean historical data may be incorporated into the input, 
but such data are neither stored nor generated internally. 
Thé specific assumptions and constraints used in each sub- 
program are described in the subsequent chapter which deals 
Eth that sub-program. 

General worksheets for PEAPS have been designed to 
guide the operator through the sequence of steps required 
to generate output (Appendix B). 

Operator note sheets also have been written in 
outline form to describe the use of the individual sub-programs 
of the system and to assist the operator when problems and 


questions arise (Appendix C). 
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ТТ. SPEED OF SOUND PROFILE (SSP) GENERATOR 


А. PURPOSE 

The speed-of-sound generator produces a plot of sound 
speed versus cepth for inputs from a bathythermograph (BT) 
trace. The data are then searched for critical points which 
are stored for subsequent use by the ray trace and transmis- 
sion loss (TL) programs. A critical point iS a point on the 
sound speed profile where the gradient of the sound speed 
changes Sign or becomes zero. All sound-speed points are 


stored for use in the ray trace sub-program. 


Be INPUTS 

Figure 1 shows the SSP portion of the operator's work-- 
sheet. The operator manually inserts the BT data from the 
worksheet into the calculator programmed for SSP/Ray Trace. 
As indicated, the data entered on the worksheet are the 
number of BT points to be entered, and the depth (ft), 
temperature (°F), and salinity (ppt) of each input point; 
Meesalinity is unknown and historical salinity data are 
unobtainable, a value of 35 parts per thousand (ppt) is 
entered. If tne bottom depth is greater than that of the 
BT trace and historical data are not available, it is recom- 
mended that the operator enter the last BT temperature as 
the bottom temperature. This will produce a profile for 
Semotent temperature from the last BT point to the bottom. 


ТЕ tables of historical salinities and temperatures-versus-depth 





PEAPS OPERATOR WORKSHEET 


(See PEAPS Operator Note Sheets) 





BT Information DTG: 2 o 
Location: N/S , E/W 
ШШ итп the calculator and plotter on. "LOAD" the SSP/Ray 


ace program cards. Put the plotting paper on the 
plotter and set the "GRAPH LIMITS" of the plotter. 


2. Press "STOP, END, CONTINUE" (this sequence is used when 
a new SSP is to be generated. 


Note 1: If you enter an erroneous value: 


а. If you have not pressed "CONTINUE", press 
"CLR X", then enter the correct value. 


b. If you have pressed "CONTINUE", you must 
re-enter all the points again; press "GO TO, 
LABEL, A", then continue with Step 3 below. 





mee Press "CONTINUE", 


ИЛ 55Р inputs: ENTER number of BT points (15 maximum) 
Press "CONTINUE" after each vaiue is entered to advance 


the program Output 


input. dut Input Speed of Sound 
Beane Depth (ft) Temp (°F) Salinity (ppt) (ft/sec) 


1 0.0 


2 


15 


FIGURE 1 
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for the oceanic area of interest are available to the 
operator, then use Section G of the Operator Note Sheets 
(Appendix G) to merge the historical and on-scene BT data 


mer improved accuracy. 


С. MODEL 

The basic physical concepts of sound transmission are 
fairly simple, but computation of the sound speed for the 
ocean is rather complex because of the varying physical 


properties of sea water. The computation of sound speed from 


| empirical relations based upon laboratory measurements began 


in the early 1950's (Del Grasso, 1952; Wilson, 1960). 
Attempts to record accurate hydrographic data while elec- 
tronically measuring the sound speed vere successful in the 
early 1960's (Hays, 1961; Mackenzie, 1961). For these latter 
experiments, a computed sound speed was obtained based on 
three measured parameters: temperature, salinity, and depth. 
The results were then compared with the electronically- 
recorded sound speeds. These observations resulted in either 
slight modification of the coefficients in existing equations 
or the development of new polynomial equations (Wilson, 1960; 
EE 1968; Frye and Pugh, 1971; Anderson, 1971; Del Grasso, 
1972). From the results obtained by these experimenters, the 


variation of sound speed versus temperature, Salinity, and 


depth can be approximately represented as indicated in Table 1. 


Approximately 99% of all sea water falls within the 


temperature range -3°C to +30°C and the salinity range 33 to 
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TABLE 1 
Variations of sound speed caused by variations in 
temperature, salinity, and depth. A variation of one 
parameter Signifies that the other parameters are held 


constant: temperature at O?C, zero depth, salinity at 


2.00 ppt. 
Temperature (°C) Sound Speed Change 
0 - 10 4.5 (meters/sec)/?C 
10 - 20 3.0 (meters/sec)/?C 
20 - 30 2.0 (meters/sec) /°C 


Salinity (ppt) 
33 - 35 1,55 (meters/see)/ppt 


35 - 37 1.25 (meters/sec)/ppt 


Depth (meters) 


0 - 10,000 0.017 (meters/soc)/meter 
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ppt (Montgomery, 1958; Frye and Pugh, 1971). Thus, 

thermal gradients are usually much more significant than 

salinity gradients in sound speed computations; so that the 

temperature-depth relation is a factor which must be 

considered in greatest detail in any sound propagation model. 
Because the Frye and Pugh equation [1] for sound speed 

is dominated by temperature terms and requires relatively 

few program steps compared with other speed of sound equations, 

it was used in the sound speed profile sub-program. The 


equation was modified to use the approximation 


for the speed of sound SS in meters/sec, the pressure P in 
bars (gage pressure), and the depth Z in meters. The 


resulting equation is 


Е 149.3 + 1.5848 х 10-25 + 1.572 х 107/322 
_ 3.46 х 1071954 + 4.5877 – 5.356 х 10727? 
-4 3 B -2 2 
о Toner + 1.19(S=35) + 9.6 x 10 (8-35) 


+ 1.354 х 103122 - 7.19 х 10-9727 


0 e (s - 35)7 


where SS sound speed in meters/sec 


Z = depth in meters 


= 
| 


= temperature in °C 


un 
|] 


salinity in ppt (parts per thousand). 
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The SSP sub-program converts the input data to MKS units 
for use in the above equation and then converts the 
resultant sound speeds to feet/sec for use in subsequent 
routines. 

The above equation is relatively inaccurate in areas 
where Salinity exceeds 36.6 ppt: Frye and Pugh claim 
accuracy only to 36.6 ppt and are presently working to extend 


the domain of their equation. 


№ OUTPUT 

The calculator prints the input data and the corresponding 
sound speeds on tape (Fig. 2), makes a plot of sound speed 
versus depth (Fig. 3), and refers the operator to Operator 


Note 2 for ray trace preparation (see Appendix C). 





зар figure 2. Sample Speed of Sound 
ОИ 0ТЕ 1, Profile generator 
printer input/output. 


ы 
ЕП 
+ 


| Е б " E А І ie 35) x 
"E TEMP < TEM 
> ! : үр 


SALIH SAL TH 1 EIN 


E 













1 
a 
| 
1 
1 


п 
і 
Do 
“A 
4 
RR 
ai 
= 
= 
E 
E 
4. 


eek 


TEC 
= 
ТЕСЕ! 
ы 
=: 


i] 

CO 

п 

т, 

га 

+ 

(rm 
má 
u 


ri 
Lf 
кз 


та 
+ 


25 Ans 358, Eta 35,.080x 


ES — — nan | SPD 


4941.78 4898,36 4077.91 
1 


ЈЕРТН 


Ізі 


H | | DEPTH 


. = . 
~ ~ - - == - 


Tan 
jn 
(a 
а 
E 
pon 
— 
{Б 
L 
1 
"fa 
E 
= 
ut 
Zt 
— 
ке 
т 
TES 
1 
1 
' 


іа! 
т 
en, 
一 
а 
эа 
— 
= 
С 
ду 


ГЕНЕ 


== 
in 
EE 
— 


ШЕШЕ 


а ЕН dr, Эйе 45. Hk 


SHLiH ШЕ 


eo 
1 
f 
== 
— 
= 
et 


HLIH 


A a в ІЗ А Bb a = 5 || % ae "a E р E E p 
^ E il - F 1 о [^ Т; 
4 9 4 А z Е, 4 d = = e ш = 4 = x Е; . г; 1 


ПЕРТН DEPTH | DEPTH 


ES 35 AAs ~ = 


<. 
(n 
t 
BIG 
iit 


Ж 


Eos 
(fü 
ара 
ЖИ 
—-l 
su 
E 
DT. 


ПЕРТН DEPTH HEP TH 


! 
| 
E 
P 
+ 
1 Li 


ioe 
=> 
— 
#..# 
=> 
u 
= 
=== 
a= 
-—— 
t 
8 
Re 
: 
ЕП 
— 
Bur 
1 
be 
ہے‎ 
د‎ 
таш 
ME 


ГЕМР | 2" ТЕМЕ ^ [4 TEMP 


ta'i 
ae 
п 
= 
ps 
КЕС 
— 
1 
4. 
г. 
po. 
E 
= 
1 
п 
La 
, 
X 


SALIH SALIH SALIH 


FEE OPT SET 


Er 
күз 
E. 
n 
co 
- Т 
t 
п 
= 
La 
DR 


1 
1 


ta 
р EL 
іл 
и 
La 
гі ез 
DOR 
4 
س‎ 


DEPTH DEPTH DEPTH 
| | E jü?5,HB8* Уса, й 


b TEE © TEMP 


m 


(£n 
zu 
a 
ma 
1m 
+ 
ES 
ПАК 
а 
ЕП 
Emi 
+ 
% 
~ 
21 
a 
E 
— 
1 
+ 


SALIH SALIH‏ ك 
Е e So MELIUS cic НІЛ‏ 


Mir el, dorr ая 


23 





*¿ndano 133901d дозедацэЗ 
ӘТТІ924 рипос JO pasds aTdues “€ exnStj 


CAZ | 


~ 


“ап Eh 


SII .L.4] Cc 


HH HH HH HH 


(az. І 


Se) 


“За І 


[ZAl X L4] Hic3Q0 


24 





III. RAY TRACE GENERATOR 





А. PURPOSE 


A ray trace sub-program included in PEAPS makes the 






following tests and decisions: 

]. inspect the speed-of-sound profile to determine 
whether the characteristics necessary for convergence zone 
propagation are present or not; 

Zee af CZ conditions exist, determine the maximum and 
minimum CZ ranges, where minimum CZ range is defined as the 
range from the source to the inner limit of the CZ annulus 
Р maximum CZ range is defined as the range to the outer 


limit of the CZ annulus (Fig. 4), and in any event; 


3. determine the minimum bottom bounce angle. 


BR MODES 

The ray trace sub-program has three modes of operation 
designated by 0, 1, апа 2. 

Mode 0 automatically determines and plots the limiting 
rays of the convergence-zone ray-family and stores the 
corresponding ranges. 

Mode 1 is the same as Mode 0, but omits the ray trace 
piot. 

Mode 2 gives a plot of the operator's input ray-family 


followed by a Mode 1 run. 


2 
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EL INPUTS 

The operator enters the Mode number and source depth 
from the ray trace portion of the worksheet. If Mode 2 
is selected, he must also enter twice the maximum and 
minimum emission angles and twice the desired increment in 


emission angle between successive rays! (ETO ST. 


I MODEL 

The ray trace generator uses Snell's Law in a small 
angle approximation. 

Snell's Law for rays can be written as 


Ат В (511092 - $1101) 


AZ 


В(со$91 - cos02) 





where the geometry is suggested in Fig. 6. Elimination of 
02 gives the range-depth equation in terms of the initial 


angle of depression, 


1. 
ED E RIL = (sinol + 55,413 (1) 


| ice Section B of Note 2 in the Operator Note Sheets 
| Аррепа1х C). 









| >: change the plotting paper on the plotter. Enter 
the ray trace program MODE indicator. 





Ў = Automatic compute and plot 
MODE: Ше xabrcocompute, no plot 
2 = Manual inputs for plotting, Automatic 


compute. 


If a ray trace input error is made, press "GOTO, 
BABEL, B, CONTINUE", then input the MODEL indicator 
and ray trace data again. (This sequence is used 
cOMgenerate a mew trace using the old SSP data) 


Press "CONTINUE" 


ША Ray trace inputs: 
Press "CONTINUE" after each input to advance the program. 


Source Depth: | Crt) 


If Mode 2 is used: 


Minimum Ray Angle x 2: (degrees) 
Angle Increment x 2: (degrees) 
Maximum Ray Angle x 2: (degrees) 


| 3: If, while in Mode Й and after the bottom limiting 
ray has been drawn, a channelled ray is drawn to 
the plot's range limit; see Operator Note sheet 
section Note 2E. 





Change the plotting paper on the plotter. Press 
"STOP, END". LOAD the TLOSS/PD program cards. 
Press "STOP,END". (The "STOP, END" sequence is 


used when a TLOSS/PD input error 15 made, or a new 
TLOSS/PD plot is to be generated based on the p 
SSP/Ray trace data.) 


5. TLOSS/PD inputs: 
Press "CONTINUE" after each input to advance the program. 
Press "CONTINUE" 


Source Frequency: (Hz) 


Apure of Merit (Ром): (ав) 

System Sigma: (use 10 if unknown) 
Me eight: р 1 23456789 10 11 12 (ft) 
Receiver Depth: fe) 


Bottom Type: 12 3456789 


FIGURE 5 
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Illustration of ray curvature within 
a water column. 





sxpanding AZl as a MacLaurin series about Ar = 0 gives 


1 - (Ar)? M $51101 BU 
AR(cos01) 


2R^ (cos01)? 


А72 = +ап01 (Ах) + 3 


1 ~ 4 (51061) 2 


3 (Ar) * паса x (2) 
8R^ (cos01) 


+ 


Substitution of the trigonometric series representations 
of tan x and cos x into the first two terms of equation (2) 


gives 


(1 + 1.5612) 


2 
5H (Nr) > < (3) 


912 
623 = 01(1 + 4) (ar) + 


а 


To conserve program space and reduce the computation 


К PEAPS uses a Simplified form of Eq. (3). Use of the 


small angle approximations 
| 


| 
CD 


sind 
ага 


cos6 


|| 

е 

! 
N 


where 9 is in radians, yields the simple approximations 


AZ 


| 2 
(81) (ar) + AL (4) 


and 


A0 


Ar/R . 
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Ес. 14) to accurately represent Eq. (3), restrictions 


mist be imposed on Ar and 61. In order that omitting а 


ntroduces no more than 103 error, 01 must be less than 
|. 33 radians. To determine the limits on Ar, examine the 
ТӘЙЕ ОҒ "worst possible" conditions: 

1. А maximum gradient of 0.5 ses + and a minimum sound 
speed of 4700 ft/sec are assumed so that the minimum radius 
£ curvature is 9400 ft. 

2. The maximum allowable angle of depression of 0.35 
:адјап5з (209) is taken. 


^o ensure that no layer is missed or poorly represented in 
іне ray trace, five depth increments are required within any 
| 

layer Of depth D. Then, the first-order term in Eq. (4) 
Las the equation 


| 
| D _ 
| E < Sel ^ OF в (5) 


| 


iS an approximate bound on the allowed Ar within a layer 
| 
2Ё depth D. 
лав Еће truncation of Eq. (2) to second order in Ar 


| 
introduces no greeter error, the third term of the equation 


| 


nust be much less than the second, 








(Ar)? «« —— (Ar) ^‏ 1هن 
AR" (cos601) AR(cos61)‏ 
\ 
R(cos81)*‏ 
En << == БВ‏ 
== $1001 
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ГЕСАТ 15 required to be at least one order of magnitude 


smaller than 2.6R, then 


Ar < 0.26R < 2444 ft. 


However, Since for small 61 the third term becomes 


arbitrarily small, it is also necessary to require that the 


fourth term of Eq. (2) be much less than the second, 


| 1% 
1 u. eye eo 1 : (Ar)? | 
8R~ (cosé@l) 2R(cos81l) 
or 
2 DT 
ШЕ << (2ғ) (сов0Ө!) [1] + 4(ејлп01) ]“ = 218. 


If Ar is to be at least one order of magnitude less 


Uan 2.1R, then 
ЕС 0-21К < 2012 ft. 


To estimate the error in calculating AZ from Eq. (4), 


Î subtract Eq. (4) from Eq. (3): 


3 2 
23 - AZ = 2 Beer)? (6) 





4R 


(АЈ) 
М 





The fractional error is given Бу: 


423 - AZ _ 0.014 + (9.56 х 106) (дк) 
Bais E Tc 10 (Ar) 


Away from a turning point, for large angles where the 

E rence defined by Eq. (6) is large, the error is 
approximately 4% for small Ar and approximately 8% for large 
Bf the error is required to be no greater than 5%, then 
Ar must be no greater than about 1000 feet (near a turning 
point, AZ3 - AZ approaches zero very rapidly, so that errors 
in depth are negligible for small angles of elevation or 
depression). 


Based on these error requirements, the simple equation 
L 
ee = 5(D) ° 


was used for Ar and D in feet. For the ranges of D employed 
in PEAPS the restrictions on Ar set forth above are not 


Seriously violated (Table II). 


ES OUTPUT 

ІІ convergence zone (CZ) conditions are not present, 
"NO CZ" appears on the printer. In all cases "SEE NOTE 3" 
and "END" appear, indicating the end of the SSP/Ray Trace 
‚programming sequence, and referring the user to Operator 


Note 3 (Appendix C). 
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TABLE II 


Comparison of PEAPS and theoretical range increment values 





D IDO, 10009) Ar = 5 [р] 
50 2805 а 
100 57.0 50.0 
| 200 114.0 70.7 
| 400 228.0 100.0 
| 800 456.0 141.4 
| 1600 912.0 200.0 
| 3200 1000.0 282.8 


6400 1000.0 400.0 








Graphic outputs depend upon the Mode number chosen. 
Figure 7 is an illustration of a Mode 0 гоп when CZ conditions 
are met. Figure 8 shows a Mode 2 run. These pictorial 
Jutputs may serve as training aids in either a classroom 
»r a shipboard environment to demonstrate the dependence of 
ray paths on the vertical sound speed structure. Mode 1 


las Mo graphic ray trace Output. 
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IV. TRANSMISSION LOSS GENERATOR 


ATH SELECTION 

Transmission loss is computed for paths which meet 
the following criteria: 

Surface Duct Path (SDP) 

a. The sound speed profile (SSP) contains a sonic 
layer. A sonic layer is defined as the layer above the 
first negative sound-speed gradient (Fig. 9). 

b. The source and receiver are both located within 
the sonic layer. 

C. Sound remains trapped in the sonic layer. This 
| trapping fails to occur at very low frequencies, where the 
wavelength has become too large to "fit" in the duct. Using 

a sound speed of 5170 ft/sec and a Sound speed gradient of 


0.017 = l in the sonic layer, we have from Urick [2] 


Oe Minimum duct thickness (H in? in feet which will 
trap sound with wavelength A in feet. Urick states that 
| this expression does not represent a sharp cutoff. To 

| reduce prograiming steps for the calculator, the dn 


, equation was simplified to 
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CD1 = Surface 
CD2 = Sonic Layer Depth (SID) 


figure 9. Illustrations of Sonic Layer Depths (SLID). 





Табје 111 compares the frequencies and corresponding minimum 
duct thickness as determined by PEAPS, Urick, and ASRAP ТТТ 
@eemectcion tables, 
2. Reliable Acoustic Path (RAP) 
RAP conditions always exist except when the source 
and receiver are both located in a sonic layer. 
3. Bottom Bounce Path (BBP) 
Bottom Bounce Path always exists. 
Convergence Zone Path (CZP) 
CZP is computed if an SSP with a sound speed excess 
is present (Fig. 4). 
5. Deep Sound Channel Path (DSCP) 
The Deep Sound Channel Path exists when the source 
and receiver are both located in the deep sound channel and 


the sound speed profile can support CZ propagation. 


Bae UTS 
The operator loads the TL/PD program and enters (from 
the TL/PD portion of the worksheet) (Fig. 5) source frequency 
(Hz), Figure of Merit (GB), System Sigma (dB), wave height 
(ft), and MGS bottom type. Although the figure of merit 
ІЗ Ж Сеп sigma are used only in the PD generator, they 
ire entered at this time to simplify operator procedures. 
The transmission loss sub-program uses the following 


ralues from the ray trace generator: a CZ indicator, the 





iaximum and minimum CZ ranges where CZ conditions exist, 


зе minimum bottom bounce angle, and the depths and speeds 
| 
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T ASIDE TII 


Comparison of minimum duct thickness 


ов екше — Hen CEE) 





Frequency 
6) ИТЕ) PEAPS ОСК ASRAPS III 
SE З L89 7 776.2 750 
104 49.7 дО. 495.3 475 
208 24.9 SUI 7 312.4 300 
588 8.8 146.6 155.2 150 
1080 4.8 95-9 104.3 100 
| 3085 7 46.4 5292 50 








of sound at the critical points of the sound speed profiie 


(sonic layer, deep sound channel axis, etc.). 


С. TRANSMISSION LOSS MODELS 


The program looks at each patn selected and determines 


transmission losses according to certain simple models. 


MS used in this section are those used in the FORTRAN 


version of the PEAPS program (Appendix D). 


ШІ. 


where R 


Surface Duct Path (SDP) 


Eure 10 illustrates SDP. 


TL 


|| 


spp = 20log(R) + AR + DFL*R for R< RT, R< RS 
ЕЕ. - 101098) + AR + DFL*R + 10log(RT) 
for ЕТ < R< RS 
SDP — ВЕЛ tank + DBL*RE JOlog(RT) + SSL*R 
LOrE Е, К > ЕТ 
spp "^ 201o0g(R) t AR * DFL*R 4 SSL*R for RS < R < RT 


шае horizontal distance from the source to the 
receiver in yards 


Бие epbsorption coetficient [3], 


mie ^ 49r^* 
l BE Е 2 
NE^ — 43004 


=> 








UO (ав/уд), 


frequency in kilohertz (KHz); 
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DFL = diffractive leakage out of the duct based on 
Normai Mode Theory [4 and 5], given by 


744p 1: 66570.33(0p52)7? ж 1073 (ав/уа), and 
ы ЫЕ = CSa c 
G = the below-layer gradient = CB3 — GD2 ! 


CD2 = the sonic layer depth (ft); 


CD3 = the depth below CD2 at which the next positive 
gradient begins or the bottom depth (ft); 


CS2, CS3 = the sound speeds (ft/sec) which correspond 
to CD2 and CD3, respectively 


RT = the transition range (the range at which the duct 
is considered full and cylindrical spreading begins) 
; _ Сри 
in yards [2] 7 бъуссов (59757) 

SV = the speed of sound at the source (ft/sec). 

RS = the skip distance (the renge between consecutive 


Surface reflections of the deepest trapped ray) 


2arccos (CSi/CS2) 


in yards [2] = “ТАРА” (yd) 


re the inverse of the radius of curvature of a 


ray within the surface duct = CD/SD. 
SD = the depth of the source of sound (ft). 
GD = the effective gradient in the scenic layer 

ШЕ - csl 
CD2 = CDi 
0 apo - 71 

SSL = sea surface loss = EE (ав/уа) 
В = РАН 
Н = wave height te) 


The expression for SSL was developed from a power regression 
analysis of data presented in Reference 6, and is valid for 


values Of B between 0 and 3. 





2. Reliable Acoustic Path (RAP) 


prom Fig. 11, 


TLRAP аа шэсць ДК, 
where K is the slant range from the source below the layer 
to the Surface at range R, in yards. 

The maximum RAP range = RL + RS/2, where RL is the 
cutoff range, defined as the horizontal distance from the 
source to the point at which a ray just grazes the bottom 
Of the sonic layer (or surface if no sonic layer exists), 
and RS is the skip distance (described previously in the 
Surface duct model). 

DEA Bottom Bounce Path (BBP) 

Four separate bottom-bounce paths are possible between 
the source and receiver: two with one surface reflection, 
one with two surface reflections, and one with no surface 
reflections (Fig. 12). Incoherent superposition of these 
rays at the receiver results in a 6 dB reduction in the 
transmission loss, For emission angles greater than about 
15°, the radius of curvature of a ray is usually large enough 
to assume straight-line propagation of the rays from the 
SOurce to the bottom and back to the surface. A ray from the 
SOurce will strike the bottom at an angle of incidence 


estimated by 
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5 I CDG = SD 
THBL = tan EH 
where 
__ 8 (Срб ~ SD) 
ХВ = N ODE Ер) SD (yds) 


is the horizontal distance between the source and the range 
of the first bottom bounce, CD6 is the bottom depth (ft), 
and N is the number of bottom reflections. 

The distance traveled by the four rays which strike 


the bottom once can all be approximated by 


R 
cosumiBr); (595) 


for reasonably small source and receiver depths or large 
range. Therefore, an approximate transmission loss can be 


тел а5 


К BER 


EE о E E A DURS Zo MOL NK 
TL cl cette.) ° cos (THBL) 


* Б 
вв T N BL 6 


where the subtraction of 6 takes into account the 6 dB gain 
caused by multiple bottom bounce paths, and BL is a bottom 
loss term which was determined through regression analysis 

of bottom loss data presented in Reference 7. Specific values 
for BL are determined by a sorting routine which uses the 
input frequency and bottom type.  THBL is the angle at which 


шатаў strikes the bottom. 
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4. Convergence Zone Path (CZP) 


7 


See Fig. 4. 


= + 
Thay РОВ) t AFR 


when R is outside the CZ annulus and 


ШЕ: - 201о4(8) + AYR - GCZ 


when R is bounded by CZSTRT and CZSTOP. The CZ annulus is 
defined as the width of the zone, CZ2STOP - CZSTRT, in yards. 
GCZ is the reduction in transmission loss caused ру the 


increase in sound pressure level (SPL) which results from ray 


focusing within the ray family forming the CZ annulus. 


ОСА = BP ехр(-вр2/2СВ) 





ЕВ 
_ ВАМСЕ 5 
о зах (GB! 
ав = 0-36 (47?) 
GCZMAX 


RANGE = R - CZMIN (yds) 


Х = CZ annular width (yds) 


4 


GCZMAX = 15.4545 - 4,545 * 10 “Х (yds) 


The equations for GCZ and GCZMAX are empirical and were 


derived from analysis of numerous transmission loss curves 


У 


obtained as output from the FACT model [5]. The analysis 


> 
мо 





established that the gain distribution within a convergence 
zone closely resembled a Rayleigh distribution and that the 
gain tended to a maximum at approximately 18% of the annular 
width. Thus, GCZ was modelled as a Rayleigh distribution 
with the peak Occurring at a distance from CZSTRT equal to 
ШОО the annular width (Fig. 13). 

5. Deep Sound Channel Path (DSCP) 

See Fig. 14. 


TL = lOlog(R) + дхв + 10104 (ЕТЕ) 


DSCP 


where RIP is the deep sound channel transition range, 


Ze CZMIN ,D. 
RTP = Б (5) 
р = CD3 - CD2 in feet, 
CD3 = depth of the sound channel axis (ft), 
7 = SD - CD2 (ft), 


CZMIN is in vards. 


The transmission loss term is valid only for R » RTP, 


although it is used for all R when this path is present. 


D. COMBINED TRANSMISSION LOSS 

Since transmission loss caused by bottom bounce is always 
Present, the total transmission loss (STIL) iS initialized 
with TLBR after each increment of range. Transmission loss 


for subsequent propagation paths meeting their respective 





figure 13. Convergence Zone Gain Model 


GCZ 


2 


ШОТ Ото - CZSTRT) 


GCZMAX = 15.5 - (0.00045(CZSTOP - CZSTRT}) 
GB = 0,36 / GCZMAX^ 
RANGE = pb - CZSTRT 


RP = (RANGE x VGB) | 
(O. 18(CZSTOP - CZSTRT)) 
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existence criteria are computed for the Same range and 
combined with STL as follows: If 0 signifies the combining 
of decibels (dB) for independent (phase incoherent) signals, 
then 


STL O TL(i) = STL = A 


where TL(i) is the transmission loss for the specific 


Propagation path (CZ, SDP, DSC, RAP, as they exist). 


UCET 


EE Ооо +10 77) 
\ 
5 - Баш E маў. 
E. РОТРОТ 


The program plots the combined transmission loss (STL) 
in dB versus range in nautical miles (NM), with a reference 


ШИЕ Т От Тһе іпрпіс figure of merit (Fig. 15). 


INS HORTCOMINGS OF THE MODELS 

Since the models were developed for use in a relatively 
small programmable desk-top calculator, limitations in the 
available storage and the number of program steps have 
resulted in certain shortcomings: 

1. Phase coherence is neglected. 

2. Sea surface scattering is not taken into account in 
calculating TL for the convergence zone model. 

3. If the source is in the sonic layer and the duct 
thickness is less than Hain’ propagation in the layer is 


neglected. 
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Д. The use of straight-line propagation will induce 
some error in the angles computed for bottom bounce 
Parameters and in the transmission loss computation. 

5. Convergence zones are considered when any sound 
Speed excess iS present. 

6. The convergence zone gain term, GCZ, is found from 
an empirical formula: there is no physical justification 
for the use of the Rayléigh distribution. 

7. The model for the deep sound-channel overestimates 
transmission loss for ranges less than the (deep sound- 
channel) transition range. 

8. For reliable acoustic path, steep-gradient curvature 
is not considered; transmission loss is computed for straight 
line distances between the source and receiver. 

9. The absorption loss term uses Thorp's formula for 
all frequencies. As shown in Fig. 16, the formula is good 
tor frequencies greater than 20 Hertz. 

10. The diffraction loss term, DFL, is the same term used 
in FACT [5]. “The following is a qutte from Reference 5: 

A preliminary analysis of this model (Labianca, 1972) 

has shown that the leakage term does not correspond to 

the mechanism which dominates duct leakage and, in 
particular, has the wrong dependence on the below-layer 

(thermocline) gradient. An intensive investigation of 

all available surface duct models has been completed 

ВЕРСТУ NUC (Naval Undersea Center) for LRAPP (Long 

Ande Acoustic Propagation Project) (Morris, 1974). This 

Study compared several models with each other and the 

SUDS (Surface Duct Sonar) transmission loss data 

(Cummins, 1972). Since none of the candidate models 

emerged as clearly superior, AESD (Acoustic Environmental 


Support Detachment) is currently studying possible 
modifications to the surface-duct equations. [5] 
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FREQUENCY IN KHZ 
Measured data taken from the North Atlantic (Thorp end 
Bernier, 1959 and 1962; Sussman, MacDonald, and Kànahis, 
1963), Mediterranean (lallement and Waterman, 1963; Leroy, 


1964), and Northern Pacific (Greer and Polam, 1954; 
Sheehy and Halley, 1957) [3] 
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11. Convergence zone parameters are determined et the 
surface. Errors exist when source and/or receiver are at 
appreciable depths. 

12. When CZ conditions are not present, the model 
assumes there is no deep sound channel. 

13. Leakage of sound into a shadow zone is not con- 
sidered. If thereare no CZ conditions, then the only 
illumination of the shadow zone comes from bottom bounce. 
If CZ conditions exist the shadow zone is illumineted by 
a path having spherical spreading in addition to the bottom 


bounce illumination. 
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O PROBABILITY OF DETECTION (PD) GENERATOR 


ШЕ. INPUTS 
The sub-program for probability of detection contains 
values from the cumulative distribution function (CDF) for 
a standardized normal random variable, CDF = F(x). Values 
Of x corresponding О 0775 9,55, ..., 0.95 are stored. 
The figure of merit (FOM) and system sigma are entered 
by the operator as part of the transmission loss initializa- 
tion procedure. System sigma is defined as the Square root 
of the sum of the variances of the individual terms in the 


pasSive sonar equation: 


| A 2 2 2 2 
System sigma = (02, + 07, + Op, + бур + Отт) 


where o is the deviation. 
Rule-of-thumb values for system sigma are included in 
Operator Note 3 (Appendix C) for those cases in which variances 


ле individual term are unknown. 


B. MODEL 
The detection variable X is defined as the ratio of the 
signal excess, 


= FOM = ТІ, 


to the system sigma: 


M SE 


_ System Sigma 
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The computed X is compared to the stored values cf x 


to determine probability of detection at which PD = F(X). 


E OUTPUT 


The program plots a probability of detection (PD) versus 
range directly below the TL versus range plot (Fig. 15). 
The PD routine computes PD in increments of 0.05 so that a 


smooth curve does not result. 





Е CONGHUDENG REMARKS 


Advanced development of an on-board environmental 
ction system resulted in the Integrated Command ASW 
Prediction System (ICAPS) computer program. ICAPS uses the 
FACT transmission loss program so that it requires computer 
resources comparable to those available on aircraft carriers. 
This development effort is continuing with emphasis upon the 
development of a smaller version of FACT. A viable alternative 
would be an augmented version of PEAPS. 

Fleet introduction of an environmental prediction system 
need not be delayed until a foolproof system is available. 
With PEAPS available, even in its present limited form, 
operational evaluation of an interim system with some known 
shortcomings is entirely feasible. A continuing dialogue, 
carried on between the developmental agency and the operating 
forces, would facilitate implementation of system improvements 
based upon fleet experience. 

The operating forces have never had a tactical environ- 
mental prediction capability available at the individual 
platform level.  PEAPS provides an inexpensive way to evaluate 


the operational utility of such a system. 
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APPENDIX A 


PEAPS HARDWARE DESCRIPTION 


For this program the Hewlett-Packard Model 9810A 
calculator must be configured with the following options: 

Option 001: Total of 111 data-storage registers; 
Option 003: Total of 2036 program steps; 
Option 004: Printer. 

Read-only-memory units (ROM's) required are: 
Mathematics ROM - hp - 11210А 
Printer Alpha ROM - hp-11211A 
Plotter ROM - hp = 11215A. 

The Hewlett-Packard Model 9862A Plotter is required to 


produce the graphical output. 
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APPENDIX B 





Ме following pages are samples of the Operator Worksheet 


eloped for use with the PEAPS program. 
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ВЕРЬ ОБЕБАТОВ WORKSHEET 


f 


(See PEAPS Overator Note Sheets) 
 Гптоггасіоп | DTG: 7 
оса јоп: VS ae E/W 


EN "n the calculator and plotter on. "LOAD" the SSP/Ray 
Кассе Program cards, Put the plotting paper on the 
EN uter and set the "GRAPH LIMITS" of the plotter. 


Beeeecess "STOP, END, CONTINUE" (this sequence is used when 
ES SSP is to be generated. 


ФЕ Tf you enter an erroneous value: 


ЕЕ ШІ уосп have not pressed "CONTINUE", press 


-e m 


"eGR xX", then enter the correct value. 
feat you have pressed "CONTINUE", vou must 

me=senter all tre points again; press "GO TO, 

Mea) A, then continue with Step 3 below. 


сә 


BEEN ess "CONTINUE". 


ОН: ENTER number of BT points (15 maximun) 
NS CONT INCE” after each value is entered to advance 
Be. program 





: OUTO Е 
при прат Inpa. > DESARRO Sounc 
ME Depth (ft) Temp (°F) Salinity (ppt) (Et/ SC 








E N NA 

2 - س‎ 
10 E 
15 
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шаша: change the plotting paper on the plotter. Enter 
the ray trace program MODE indicator. 


f = Automatic compute and plot 
ODE: о: < CONPULe, no plot 
TTI INPDULS for plotting, Automatic 


compute. 


ШҰҒа гау trace input error is made, press "GOTO, 
LABEL, B, CONTINUE", then input the MODEL indicator 
and ray trace data again. (This sequence is used 
Ше аспегасе a new trace using the old SSP data) 


Press "CONTINUE" 


Sy 


Пау trace inputs: \ 
Press "CONTINUE" after each input to advance the program. 


Source Depth: (MEE) 


косе 2 is used: 


Minimum Ray Angle x 2: (degrees) 
Angle Increment x 2: (degrees) 
Hum Ray Angle x 2: (degrees) 


I ТЕ, while in Mode Ø and after the bottom limiting 
ray has been drawn, a channelled ray is drawn to 
the plot's range limit; see Operator Note sheet 
Bection Note 2E. 
шшаша the plotting paper on the plotter. Press 
IE МО". LOAD the TLOSS/PD program cards. 
Press "STOP,END". (The "STOP, END" sequence is 
E cn a TLOSS/PD input error is made, or a new 
ШЕРЕС/Рр рос лв to be generated based on the old 
SSP/Ray trace data.) 


ИШОСС/РГ іпріисс: 
EN CONTINUE" after each input to advance the program. 


Es "CONTINUE" 


Source Frequency: \Н2) 

Figure Of Merit (FOM): (ав) 

System Sigma: (исе 10 if unknown) 
Сане: 91234656789 1.0 11 12 (ft) 
Receiver Depth: (EE) 


om Туре; 123456789 
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If a new ray trace is to be generated using the 
old SVP data after the TLOSS/PD program has been 
run, Press "STOP, END". LOAD the SVP/Ray trace 


готам cards. 


Note 4: 


Press “GOTO, LABEL, B, CONTINUE", then see Note 2 
above. 
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APPENDIX C 


Pepe oe OPERATOR NOTE SHEETS FOR THE HP9810A 


SSP/RAY TRACE PROGRAM 

The SSP/RAY Trace program 1s recorded on two short 
magnetic cards, each card having two sides. The program 
is inserted into the calculator as follows (it is assumed 
the calculator is turned on and paper has been placed on 
the plotter). 

1. Press: END, LOAD. 

2. Insert the cards starting with side 1 Sequentially 
through side 4. 

3. Press: END, CONTINUE, 
After the SSP axes have been drawn the printer will meh cate 
NOTE 1", referring the operator ot Note 1. 

Note 1 

A. The sound speed profile (SSP) generator will accept 
Memeo fifteen (15) points. These points are entered from 
the PEAPS Computer Operation Worksheet starting with the 
first point being the surface (depth = 0) temperature and 
Emily. The data are entered as follows: 

| 1. Press: CONTINUE | 
itive ter will request "NO. PTS". Enter the 

number of BT points to be entered as indicated on the 


worksheet, 
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2. Press: CONTINUE. 

The printer will request DEPTH. Enter the depth 
feet of the point for ШЕН data are being entered (enter 
or point 1). 

Express: CONTINUE. 

The printer will print the depth entered and 
request TEMP. Enter the temperature in °F for the point 
Of interest. 

4. Press: CONTINUE. 

The printer will print the temperature entered 
and request SALIN. Enter the salinity in parts per thousand 
for the point of interest. 

5. Press: CONTINUE. 

The printer will indicate the salinity entered 
and the speed of sound in feet per second ccrresponding to 
the depth, temperature, and salinity values immediately 
preceding the sound speed value. The plotter will plot the 
speed of sound indicated on the printer at the corresponding 
depth. If the data entered are from a point less than the 
total number of points to be entered, the printer will indi- 
cate DEPTH. Enter the depth of the next point of interest 
and continue from step 3 above. If the date entered is for 
e last point, the program will discontinue plotting and 
memeranter will indicate "SEE NOTE 2". This refers to 
Operator Note 2 below. 


Dt the BT trace does not reach the bottom: 
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l.. Use at least two points from historical data 
and enter these as the last two points. The more points 
entered, the more accurate the plot if the bottom is con- 
siderably deeper than the BT trace depth. See Section G 
below for merging data. 

2. If historical data are not available, enter the 
last BT temperature as the bottom temperature, the bottom 
depth as measured or approximated, and 35 ppt as salinity. 

C. If salinity values are not available from measurements 
n historical data, enter 35 ppt. 

D. If you enter an erroneous value: 

fee you have not pressed CONTINUE, press "CLR X", 
then enter the correct value. 

EN: you have pressed CONTINUE, you must reenter 
EENNMEUEEDoints again; press "GO TO, LABEL, A", then continue 
from step 2 in Section A above. 

E. If you have finished with one sound speed profile 
and want to enter a new one, press "END", then go to step 1 
in Section A above. 

F. Pressing "CONTINUE" after each input advances the 
program and enters the input point into storage. 

See tt the in-situ and historical data do not match 
closely together, merging becomes necessary. There is no 
single correct method for merging historical and in-situ 
data. Suggested procedures are to fare (draw a smooth curve) 


the two profiles together, to extrapolate the in-situ trace 


68 





down to tne historical data, or to extend the profiles in 
the area of interest until they meet. Fig. 17a shows a plot 
of in-situ SSP and historical data on the same axis. It is 
 есгей that this be done for all SSP plots if historical 
data are available. A trace of the historical data need 
not start at the surface. A depth 500 feet above the last 
BT depth should suffice. Figures 17b and 17c show examples 
of extending the two traces to achieve a merge. In any 
event, a smooth blending of the traces should be attempted. 
In some instances it may be necessary to try more than one 
technique in order to achieve a reasonably correct facsimile 
of the sound speed structure. 

Pate 2 

This Note deals with procedures for the ray trace 
generator. 

A. Change the plotting paper. ; 

B. Enter the ray trace program Mode indicator (0, 1, 
E 2). 

Mode 0 gives an automatic search for and plotting 
Of convergence zone critical rays. The minimum source 
emission angle that will produce a bottom bounce, and the 
minimum (inside) and maximum (outside) convergence zone 
annulus ranges (if CZ is present) are stored. 
Mode l gives the same computation as Mode 0 except 


ће ray traces are not plotted. 
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Mode 2: the operator enters twice the minimum and 
maximum emission angles to be traced from the source along 
with twice the incremental angle between. Doubling the 
desired angles and angle increment is necessary Since the 
“сап internally halves all input angles prior to conver- 
sion to radians and storing. When the manual input plotting 
is completed, the program changes to Mode 1 to determine 
the critical parameters to be used by follow-on programs. 

С. Program Mode 0 is intended to be the primary Mode 
for fleet use. Mode 1 is intended for use when a ray trace 
is not desired. Mode 2 is intended to be used for ray trace 
plots of specific source emission angles as well as for 
determination of critical parameters for follow-on programs - 
it is foreseen as a training aid. 

D. If you have made an incorrect input and entered it 
with the CONTINUE, press: GO TO, LABEL, B, CONTINUE. 
Re-enter the ray trace input parameters starting with the 
Mode indicator. 

Е. it is possible, under certain sound velocity profile 
conditions, for the ray trace program in the automatic modes 
to trace a ducted ray to the end of the plot while searching 
for the convergence zone critical ranges. If this situation 
should occur you must check the stored convergence zone 
ranges as follows: 


ЕЕ Х<( ), 008, +, 2000, +. 
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This calls out the minimum convergence zone 
menage and converts it from yards to nautical miles. 

2. Compare the value in the Y-register with the 
minimum convergence zone range indicated on the ray trace 
mot. If they differ considerably, 

a. Enter the minimum convergence zone range 
Beem the ray trace plot; 
b. Press: +, 2000, x (multiply), Y>( ), 008. 
This converts the ray trace plot minimum 
CZ range from nautical miles to yards and stores it in the 
appropriate storage register. 

Press Xel ), 009, +, 2000, +. 

This calls out the maximum CZ range and converts 
it to nautical miles. 

4. Compare the value in the Y-register with the 
maximum CZ range off the ray trace plot. If they differ 
considerably; 

a. Enter the maximum CZ range from the ray 
Brace plot. 
Би Ргевв: +, 2000, X (multiply), Y-( ), 009. 
This converts the ray trace plot maximum 
CZ range from nautical miles to yards and stores it in the 
appropriate Srorage register. 
5. Proceed with the transmission loss/probability 


of detection program. 
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Е. ТЕ you want to do another ray trace and plot using 
the same na speed profile but different input parameters: 
1. Press: GO TO, LABEL, B, CONTINUE. 
2. Enter the new ray trace inputs starting with the 
Mode indicator. 

КОО oes NOTE 3" and "END" will be printed on the output 
tape and the calculator's X, Y, and Z registers will illuminate 
when the ray trace portion of the program is completed. 

fee Oo CZ, "SEE NOTE 3", and "END" will be printed on 
the output tape to indicate that the program has ended and 
there 1S no convergence zone present. 

ЕЕЕ NO CZ", "SEE NOTE 3", and "END" is printed out 
and you would like a ray trace plot, you will have to use 
Mode 2 and manual emission angle parameters to get a plot 
(see E above). At the end of the Mode 2 plots "NO C2", 

"SEE NOTE 3", and "END" will be printed on the output tape. 

J. If another sound velocity profile is to be entered 
for subsequent ray trace plots, press: END, CONTINUE, to 


start the entire sequence again. 


T3 





TRANSMISSION LOSS/PROBABILITY OF DETECTION PROGRAM 


Note 3 

A. Change the plotting paper. 

В. The TL/PD magnetic cards are loaded in the same 
Manner as the SSP/Ray Trace cards. 

С. Press: STOP, END to initialize the program function. 
This sequence is used when a TLOSS/PD entering error is made, 
or a new TLOSS/PD plot is to be generated based on the same 
SSP/Ray Trace generated data. 

D. Source frequency is to be entered in Hertz (Hz). 

E. Figure of Merit: ensure that you use all the same 
E ot reference (uPa, or ub). DO NOT add 100 to the 
computed FOM. 

F. System Sigma 

If the variances of the sonar equation teris are 
known, enter the square root of the sum of the variances. 

If the variances are unknown, choose one of the following: 

l. Use 6 dB if ambient noise measurements have been 
made, and submarine speed and type are known. 

2. Use 8 dB if ambient noise is estimated from a 
forecast, submarine speed is known to within 3 kts, and 
submarine type is known. 

3. Use 10 dB if ambient noise is estimated from 


forecasts and submarine speed and type are uncertain. 
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С. Wave Height (IWH) 


Sea State Term Не не аме (теу 
l Calm Она 
2 Wavelets ЈЕО 2 
3 Sit 2 to 4 
4 Moderate 4 to 8 
5 Rough 8 о 13 
6 Very rough 13 to 20 
7 High 20 to 30 
8 Very high Зоо 45 
9 Phenomenal over 45 


H. A bottom type indicator (MGS) of 1 indicates little 
sound loss; a bottom type indicator of 9 indicates a large 
Bound loss. 

I. If new TLOSS/PD data is to be entered to make 
additional plots from the same ray trace data, press STOP, 
END. Change the plotting paper, then press CONTINUE. Then 
enter the new data. 

J. If a new ray trace is to be generated from the old 
SSP information after the transmission loss/probability of 
detection program has been run: press: STOP, END, LOAD. 
Load the SSP/Ray Trace program using the magnetic cards, 
change the plotting paper, press: GO TO, LABEL, B, CONTINUE, 


then follow the instructions in Note 2. 


А APPENDIX D 


Flow charts, program listings, and output from the 





310A and IBM 360 are included within this Appendix. 
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Speed of Sound Profile / Ray Trace Program 
flow charts and sample outputs 


ТЛ 


Мо 


Speed of Sound Profile / Ray 
Trace Program 


Draw the 
Speec of 
Sound 
Profile 
Axes 


LABEL A 
see Note 1 


Input the # 
of SSP pts. 


(2) 
>. 


Input Depth | 
(ft) 








Input (°F) 
Temperature 








‘Input (°/co) 
salinity 





0 
у 
2 


have 
all the SSP 
noints bee 
input? 


ЄЗ 
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\/ 


Sort out and 
store the 
critical Speed 
of Sound 
Profile points 


LABEL B 
See Note 2 


(T 


4) 


Input MODE 
indicator 









Iinatrlice 
the needed 
parameters 
for Ray 
Tracing 














| Input Source 
Depth (ft) 





MODES + 


= Auto compute and plot 

= Auto compute, No plot 

= Manual input of ray 
famiiy for plotting, 
then MODE 1 Auto 
compute with no plot 
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Sample Speed of Sound 
r Zr: Profile program input / 
printer output 
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(4) SSP critical points 
selection program 


Test = 
| Run = 
NG = 


A 


Sy = SX 
Sx = the 
next S 
Run = 






алл points 
1 st L 
"been E ted fed 


No 


Nest = 
Test + 1 








Test = ў 
store Sy 
NG = NG + i 
increment 
storage 







is 
the FLAG 
set 





Yes 
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(4) SSP critical points selection program 


Yes 






store Sy 
increment 
storage 

NG = NG + 1 
set the FLAG 






V 


store Sx 
set the FLAG 
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store Sy 





increment 
storage 


ме = ме + 1 | 


С 


Run = 
Run + 11. 












Sy = эх 


SX = the 
next S 


© 





(4) SSP critical points selection program 


О 
V 


store Sy 
increment 
storage 


NG = NG F 1 


1 


store D1 and Si 
на СРЕ and CST 


store the last 
D and S in CD6 
and CS6 


LABEL B 
see note 2 











store Sy 


increment 
storage 


NG = NG + i 
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(А) Speed of Sound Profile / Ray Trace 
Program 


is 
MODE = 0 
2 





is 
MODE = 1 
2 






Input min 
source angle 





Input angle 
increment 










Input max 
source angle 


E 


Draw Ray 
Trace axes 
(Modes 0&2 ) 





(5) 


ر 


are N, 
CZ rays (6) 
possible No. 
Б | 
Yes 
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M 
Output 
PNO CA7 | 


2 

< 

E 
"T 
UM 
^j 


N Aur) 





RANGE ENM2 





мя. 


Ray Trace axes 


(5) 


ЧЧЫ + 


ti и Ы lá 
i Ы [à 1 > 


SE. 


LES] X Jd ія 





(6) Convergence Zone Test 






Compute 
НР, 


PHPP 





Сос > CSm 


dl 6 
< 人 Yes N Е 


Compute 
SV, THP, NO 


ИР E 
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E Computation of SV, THE, HPP 





SV ,= (SD-Dm-1)X 
Зе Пеј | 


Dm-Dm-i 
Sm-1 





E 


THP = 
2SV 


“2682 


THPP = 


|2 Al 


` 056 





В Speed of Sound Profile / Ray Trace 
Program 


165 | о 








Les 


No 
THMIN=THMIN+, 59 
set plotting 
indicator 


CHMIN=TEMIN-.1° 


clear plotting 
indicator 


er >> set FLAG 


TEBB=TEMIN 
set ICZ 

set plotting 
indicator 
CZmax=R 
















Draw the ray 






THMIN = 
тир, зе 
plotting 
іпаіса- 
tor 
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(8) Ray Trace Program 







Clear 











Proven | 
Indicator МО A Surface Re- 
| | flected Ray, 
Return to Call 
3 Point 
бб = Sp ta 
E j a=) 
5 № | D(m) - D(m-1 
m= Y 
RADINV = 
Yes 


5YD(m) -D(m-1) 


e 1 © 


A Bottom Reflected 
Ray, Return to 
Call Point 
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(8) Вау Trace Program 


RDR = 
RADINV x 
DR 





Yes Xe) 


No 





Yes 


V 
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) | = to Call Point 
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Compute RMAX 
of the plot 
based on FOM 

















Input Wave 
Height (ft) 


¡Input жей 


Receiver Dpt 









Draw the 
TLOSS/PDET 
axes 





(16) 


RM 





(15) RMAX Determination 


b = FOM 






is 
80 > FOM 


? 






Yes 








No b = 80 
| 100 < POM 
Yes 
No b = 100 





| T Der 
" 





EE 
dg ° € | і. яў 


палі тоя ТЫ” 


syndur дојелодо pue «2070 Ex" pd 
(06 = WOJ) Sexe yda y ссолт әтдӛшес (от) ніала YAQA 
00 `2 Ea 
LH JAWM 
00701 
VWOIS | «zi 
02706 


JO ° SSE 2.12 


| 
| 
| 28 D 
| 8 
Wod zia 
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bas a Сы 
сма шамым 


© TLOSS / PDET Program 


is 
e «(зрэ 0р2 „#л= 
No 






“а | ما‎ = = 
IRAP 
FE uz ^ |make RL2 

<Q | (18) 
Clear | ра. 
150 у 





| set ISD | 
<i DFL = 0 „| make RS, RT, | 
аб | SSL, DFL, 
< | cloar IRAP _ | (17) 
|| No к 


CZSTRT=CZMIN | 


ICZ set | A CZSTOP=CZMAX 
2 Yes i 


_15 
1.7 
"S ч CD2 





TSD set 


clear 

IDSC NES | 2 

| set IDSC в 
make RPT и! № 


KC 
Q | (19) _ 
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(17) Computation of RS., Rf, SSL, DFL 


RS = 
CS (1; 
агссс5 с5(2 м 


асацы 8) 
D(2) - с 
(за - stay) 





Він хв] “71 
RS | 1000 














cS( 3) -cs(2)1 
CD( 3)-0D(2 | 


DFL = 744 x 


SF -5/3 x 
1600 


в71/Эхср(2)73 
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(18) Computation of RL2 


RL2 


= 
mao 


CS2(CD3--CD2 ) 


No 


RL2 = 
RL2(RD-CD2) 
4.5 


С52-С53 
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RL2 = 
RZ2(SD-CD2) 


4,5 





(19) Computation of RPT 


No 


RPT CZ iN Х | 
-/ |Сри-ср2 
SD-CD2 


8 


E 


ў 


pa 
~ 
“2 


1 TLOSS / PDET Program 


Determine the 
Bottom Loss 
factor BT 


(20) 
Determine .the 
Bottom Loss 


factors Bf, 
| Bl, ВІМА 







(21) 






Compute | 
А | 










THBB<THBL 





(22) 







store the 
Normal (0,1) 


CDF points f 
К = Е + 


DRVRK 









Бресе ТСЛЕС, 
N,» DRWRE,DRPERM 


Compute THBB, 


THBL © 
(23) 






= — 


(20) Bottom Loss Factor BT Determination 


is 
6 < BT reg] 24 8722.0 
No 
| 2| BT = 1.44 
Yes 
is A 
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(21) Determination of 
BLMA = Bottom Loss Factors 
309 BLMA, Bf, B1 
<! В1 = .35 
> > 80 = 7.12 D 
No | 
is 
! pow 29 | 
2000 < SF >— Авр = 8,21 | 
? 


No 
BLMA = 
OUR 
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(22) Computation of Absorption Coefficient (A) 


2 2 
.1(SF/1000) 40(SF/1000) 

6 + ар/уа 

& = 10007 (SF/1000)2 "4100+(SF/1000)° /у 
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(23) Preset ICZRC, N, DRWRK, DRPERM 
Compute THBB, THBL 


ICZRC 
N= 1 


DRWRK 2000 
DRPERM = 2000 












is 
ICZ set =- 
2 No 
| Yes qe 
- CS6 > SV 
Yes 
? 
У | р 
THBB = y 
THP | ЖЕ 


arccos(SV/CSS) 





THBL = 


п / 2 


ks TLOSS / PDET Program 


STL = the 


Bottom Bounce 
Path TLOSS at 
range R 









24) 





Compute the 
Surface Duct 

Path TLOSS at 
range R and 6 
LO ГІ, 








TRAP cleared 





Compute the 
Reliable 

Acoustic Path 
TLOSS at range 
R and 0 to STL 
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is 
ICZ cleared 









Yes 





Compute the 
Convergence 


at range RF and 
& to STL 





Compute the | 
Deep Sound 
Cnannel Path 
TLOSS at range 
E and © to StL 








(28) 






Plot STL 


(29) 





(24) Bottom Bounce Path TLOSS at Range i 





THBB<THBL 






ZHBL = N tan 
R(CD6-SD)/3 
WK 2CD6-SD) /3 






{ттш 
THBL х [100] 
п 








BLMA<THBLA 





Yes 


| THBLA = 


BLMA 









V 
BL = BT x N x | 


((BIXTHBLA) + 


Bg) 








STL = 20 x log 
(2/cos( THBL)) + 
A(R/cos(THEL) ) 
+ BL = 6.0 











(25) Surface Duct Path TLOSS at Range R 


TL = 
R(A + DFL) 


Yes 


No 


+ 20log(R) 





\/ 


TL = TL + 
10log(È) + 
10Log(RT) 


? Yes 





| 


TL = TL + 








Subprogram 
(see 25a) 





(25а) STL = STL 9 TL Subprogram 





AD = 10 x 


log(1*107) 


Return to 
Call Point 





(26) Reliable Acoustic Path TLOSS at Range R 


Th E 
1oc(VRD24R2 )+ | 
A(VRD 4R“ ) 


Call 
SIL = 
ӨТТ, Ф ТІ, 





Proceed to 
CZ Indicator 
Test 





No 


ED = 20 x 
ЕН )+ 
А Марса ве ) 


рат 


A AE A | 


152 













TL = 
20 log(R) + 
(A x в) 





CZSTOP « К 






DRWRK = 


10 





(CZSTOP-CZSTRT) 
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RANGE = 


R - CZSTRT 


RP = 

RANGE х Мсв 

0.18(CZSTOP 
- CZSTRT) 








6 


(27) Convergence Zone Path 










ICZRC = 
ICZRC+1 


> 


Call 
STL = 
SLL STL 


Subprogram 
(see 25a) 


TLOSS at Range R 


DRWRK 


«= 


DRPERM 


==) 


стот = 
CZSTRT + 


CZMIN 


CZSTOP = 


CZSTOP + 
С2МАХ 
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(28) Deep Sound Channel TLOSS at Range R 


TL = 
10 x log(R) 


TL + 


log(RPT) 





(A x R) 


STL e TL 


Subprogram 
(see 25a) 









TLOSS / PDET Program 









Determine the 
Probability of 
Detection at 

range R 







(30) 


Plot the 

Probability 

of Detection 

VS range .. 
(31) 





(32) 


(30) Determination of Probability of Setection 


E 


PD = 0,5 
SE = FOM-STL 


SE | 
an 





| e PD | 






vs Range 






PD = Рр-,05 
b = 5-1 


No Yes 





40Td тяаа / SSOTL әзәтӛшоо в yo eiduug (26) Т 


Eee рага 
-- 4444 


| 
| 





ш 
d 
| tu 1% 
| 
| 


(06 = HOA) не 
эзцех SA Ladd Jo 4074 atdwes (TE) + -> 2 





esuel сл Тіс JO 34oTd әташес (62) E " WL 
路 
=з 
“тя 
аы e шө ES | 


ЕН LIENS 


ааа 
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PRAPS HP 98108 
Transmission Loss/Probability of Detection 


Program Listing 
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АРРЕМОТХ Е 


PEAPS STORAGE LOCATIONS IN THE HP9810A CALCULATOR 


A. The operator enters all angles in degrees. The zero 
degrees direction for all angles is horizontal to the right. 
Clockwise angles (down) from the horizontal are to be entered 
as positive (+) angles. Counter-clockwise angles (up) from 


the horizontal are to be entered as negative (-) angles. 


B. All angles are stored and used in radians. Entered angles 
are automatically converted from degrees to radians for 


storage by the program: 


degrees = —— * radians 
180 
radians - 130 * degrees 


С. Enter all temperatures in degrees fahrenheit (°F). 


Temperatures are automatically converted to degrees 


centigrade (°C) for use in the SSP program: 


ос = (°F = 32) 


ојл 


mele 324: (= кес) 


LOS 





D. Enter depths in feet (ft). Depths are stored in feet 
and are converted to the required units as needed in the 


programs: 


meters =' 0.3048 * feet 
feet = meters/0.3048 
yards = feet/3.0 


nautical miles (NM) = yards/2000 = meters/1852. 


E. Ranges are in feet in the ray trace program and are 
automatically converted to units of yards for use in the 


transmission loss/probability of detection program. 


F. The following pages of program forms indicate the 
parameter storage locations used in the PEAPS model. The 
terms listed are defined in the Table of Symbols and 


throughout the body of the text. 
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